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Abstract 

In this paper the two-photon interactions at high energies are investigated considering different 
approaches for the QCD dynamics. In particular, we calculate the 7*7* total cross section in 
different theoretical approches and present a comparison among the predictions of the BFKL 
dynamics at leading and next-to- leading order with those from saturation physics. We analyze the 
possibility that the future linear colliders could discriminate between these different approaches. 
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The high energy hmit of the perturbative QCD is characterized by a center-of-mass which 
is much larger than the hard scales present in the problem. The simplest process where this 
limit can be studied is the high energy scattering between two heavy quark-antiquark states, 
i.e. the onium-onium scattering. For a sufficiently heavy onium state, high energy scattering 
is a perturbative process since the onium radius gives the essential scale at which the running 



coupling as is evaluated. In the dipole picture ll], the heavy quark-antiquark pair and the 
soft gluons in the limit of large number of colors Nc are viewed as a collection of color dipoles. 
In this case, the cross section can be understood as a product of the number of dipoles in one 
onium state, the number of dipoles in the other onium state and the basic cross section for 
dipole-dipole scattering due to two-gluon exchange. At leading order (LO), the cross section 
grows rapidly with the energy {a oc a^e^"^"^^-^, where {ajp — 1) = ^"^^'^ ln2 0.5 and 
Y = In s/Q^) because the LO BFKL equation predicts that the number of dipoles in the 
light cone wave function grows rapidly with the energy. Several shortcomings are present 
in this calculation. Firstly, in the leading order calculation the energy scale is arbitrary, 
which implies that the absolute value to the total cross section is therefore not predictable. 
Secondly, as is not running at LO BFKL. Finally, the power growth with energy violates 
s-channel unitarity at large rapidities. Consequently, new physical effects should modify the 
LO BFKL equation at very large s, making the resulting amplitude unitary. 

A theoretical possibility to modify this behavior in a way consistent with the unitarity is 
the idea of parton saturation Sj, where non-linear effects associated to high parton density 
are taken into account. The basic idea is that when the parton density increases (and the 
scattering amplitude tends to the unitarity limit), the linear description present in the BFKL 

equation breaks down and one enters the saturation regime, where the dynamics is described 

j 

by a nonlinear evolution equation and the parton densities saturate . The transition line 
between the linear and nonlinear regimes is characterized by the saturation scale Qsat{x), 
which is energy dependent and sets the critical transverse size for the unitarization of the 
cross sections. The saturation approach implies that the onium-onium cross section still 
grows with the energy, but the rise is slower than s°-^ Q. InRef. fl, Mueller has argued that 
these corrections become important at the values of rapidity of the order of Yj/ ~ ap-i of" 
Another possible solution, which is expected to diminishes the energy growth of the 
total cross section, is the calculation of higher order corrections to the BFKL equation. 
After an effort of ten years, the next-to-leading order (NLO) corrections were obtained 
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therein). In particular, in Ref. 
setting procedure 



6(] and the spurious singularities were solved considering renormalization-group improved 
regularizations (for a review on NLO BFKL corrections, see e.g. Ref. 3] and references 

'~\ the Brodsky-Lepage-Mackenzie (ELM) optimal scale 
9| was used to eliminate the renormalization scale ambiguity present in 
the MS-scheme calculations. The NLO BFKL predictions, as improved by the ELM scale 
setting, yields {ajp — 1) = 0.13 — 0.18. In Ref. [lOj the authors have shown that due to the 
running coupling effects the NLO corrections become important at the rapidities of the order 
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of Ynlo ^ C(s ■ Consequently, in principle, one has Yjj <^ Yj^io at parametrically small 
as- That implies that the center of mass energy at which the saturation effects become 
important is much smaller than the energy at which the NLO corrections start playing 
an important role. However, it still is an open question, since it only can be definitively 
established once next-to-leading-logarithmic contributions are fully computed for the non- 
linear evolution equations present in the saturation approaches. One possibility to check 
this assumption is the analysis of the energy dependence of the onium-onium cross section, 
which could disentangle the QCD dynamics. 

A reaction which is analogous to the process of scattering of two onia discussed above 
is the off-shell photon scattering at high energy in e~ colliders, where the photons are 
produced from the leptons beams by bremsstrahlung (For a review see, e.g., Ref. llj). In 
these two-photon reactions, the photon virtualities can be made large enough to ensure the 
applicability of the perturbative methods or can be varied in order to test the transition 
between the soft and hard regimes of the QCD dynamics. From the view of the EFKL 

nn 

approach, there are several calculations using the leading logarithmic approximation |12l . I13I | 
and considering some of the next-to-leading corrections to the total cross section 7*7* process 



13|, 



14| . In particular, the stable next-to- leading order program relying on the ELM optimal 



scale setting §1 has produced good results with OPAL and L3 data at LEP2 [l^. On the 
other hand, the successful description of all inclusive and diffractive deep inelastic data at 
the collider HERA by saturation models 15|, llGj suggests that these effects might become 
important in the energy regime probed by current colliders. It motivated the generalization 
of the saturation model to two-photon interactions at high energies performed in Ref. 
which has obtained a very good description of the data on the 77 total cross section, on 
the photon structure function at low x and on the 7*7* cross section. The formalism used 
in Ref. is based on the dipole picture (l|, with the 7*7* total cross sections being 
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described by the interaction of two color dipoles, in which the virtual 



photons fluctuate 



into (For previous analysis using the dipole picture see, e.g., Refs. [18|, ll9|). The dipole- 
dipole cross section is modeled considering phenomenological models based on saturation 
physics. A current shortcoming of the saturation approaches is that the non-linear evolution 
equations (BK and JIMWLK) were derived at leading order accuracy, resumming all powers 
of Us Ins along with all the multiple scatterings in the target, but disregarding for instance 
any running coupling corrections in their kernels (For recent efforts see, e.g. 20|). Another 
important aspect which should be emphasized is that these evolution equations contains the 
LO BFKL equation in the low density limit. 

As the NLO BFKL predictions as well as the saturation one describes quite well the LEP 
data I21I, a current open question is if the future e^e~ colliders will allow to discriminate 
between BFKL and saturation predictions. The main goal of this letter is to present a com- 
parison of the (LO and NLO) BFKL and saturation predictions for the energy dependence 
of the total cross section and verify if it is possible discriminate between these approaches 
in the kinematical range of the future linear colliders. 

Lets start presenting the main formulas to calculate the total cross section in the BFKL 
formalism, one has that the total cross section of two unpolarized gammas with virtualities 
Qi and Q2 in the LO BFKL formalism reads as follows: 
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where one has used as kinematic variables the 7*7* c.m.s. energy squared s^^ 

(p + f/)^, with "p and q being the photon momenta and the photon virtualities squared 

given by Q\ = —q^ and Q2 = —p^. Moreover, $j(z/) are the transverse and longitudinal 

polarizations gamma impact factors, which at LO accuracy are given by : 

2 / \ 2 



$t(i^) = $r(— z/) = aas 




— IV 



ri \-iv] ri i + ZZ/ 
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and 



$l(i^) = $l(-z^) = aas[ J^eJ 



r 
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The Regge scale parameter sq is proportional to a hard scale Q ~ Qi,Q2', T is the Euler 



F-f unction and e„ is the quark electric charge. 



The high-energy asymptotics of the cross sections is determined by the highest eigenvalue, 



of the BFKL equation 



3: 



0" ~ s 



ajp — 1 



At leading order it is rather large: 



12 In 2 [as/Ti) — 0.55 for as = 0.2. On the other hand, the eigenvalue of 



the NLO BFKL equation at transferred momentum squared t = in the MS-scheme reads 
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where 



XM = 2^{l) - ^(1/2 + iv) - ^{1/2 - iv) 



(5) 



is the function related with the LO eigenvalue, ip = F'/F denotes the Euler ?/^-function, the 
//-variable is conformal weight parameter Q], Nf. is the number of colors, and Qi is the 
virtuality of the reggeized gluon. Moreover, the NLO coefficient rj^s of Eq. (j4]) can be 
decomposed into a /3-dependent and a conformal (/^-independent) part [rj^{i') = r^^(z/) + 



u)] (For more details see Ref. [8|). The NLO BFKL Pomeron intercept then reads for 



cori/, 
MS * 

Nr = 3 as follows 
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(0) - -20.12 - 0.1020iVF + 0.06692/?o, 
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and 



As discussed in Refs. 



\Nf=4 



-19.99. 



23| the NLO corrections as given in the MS-scheme implies that the 



highest eigenvalue of the BFKL equation turns out to be negative and even larger than 
the LO contribution for ag > 0.157. In Ref. the authors have demonstrated that this 
situation can be improved evaluating the intercept of the BFKL Pomeron at NLO using 
the BLM scale setting within non-Abelian physical schemes, such as the momentum space 
subtraction (MOM) scheme. In this case one has 
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where 
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and 



rpconf 
^ MOM 



Nr. 



2 ' 



(10) 



with / = -2 dx\ii{x)/[x'^ - X + 1] ~ 2.3439 and ^ a gauge parameter. The /^-dependent 
part of the rMOM{i-') defines the corresponding BLM optimal scale 



Q^rJ^^M = Q^exp 
At large v one obtains 
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At z/ = one has Qblm^{0) = Q2(4exp[2(l + 2//3) - 5/3]) ~ 127. As a consequence 
of this procedure, the NLO value for the intercept of the BFKL Pomeron has a very weak 
dependence on the virtuality [8]. Moreover, it implies that (ajp — 1) = ujnlo = 0.13—0.18. 
As emphasized before, using this approach it is possible to describe the OPAL and L3 data 
at LEP2 Q. 

Another alternative to cure the highly unstable perturbative expansion of the BFKL 
kernel was proposed in Ref. 2J], who realized that the large NLO corrections emerge 
from the coUinearly enhanced physical contributions. A method, the cu-expansion, was then 
developed to resum collinear effects at all orders in a sistematic way. The resulting RG 
improved BFKL equation was consistent with renormalization group requirements through 
matching to the DGLAP limit and resummation of spurious poles. In this approach the 
kernel is positive in a much larger region which includes the experimentally accessible one. 
Based on this approach, but with a fixed coupling, Khoze, Martin, Ryskin and Stirling 

KMRS) 251 proposed recently a simpler model for the NLO corrected LO kernel (See Ref. 

251] for details). In this approach the position of the singularity (pole) is given by a implicit 
equation. 



(13) 



6 



where = Nd vros and the characteristic function x is given by 

X(z/, wnlo) = Xo{.T^) + ttsXll'^, t^NLo), (14) 
with Xo(^) = Xl defined in Eq. ([5]). The imphcit equation, eq. f|T3l) . can be rewritten as 

^NLO = tts i Xo ('^Jll - ^) + 1^ , ^^T^ ^ ^ , (15) 



, (1/2 + 07/2)2 + ^ 

where Xo* is the higher twist part of Xo, 

X'iiv) = 2i,{l) - ^(3/2 - iv) - ^(3/2 + iv), (16) 



w N 11 iVF 

AAoj) ^ ^ + 

^ ' 12 18 



(17) 



and Ai{uj) is obtained from the Melhn transform of the DGLAP sphtting function Pgg, 

'67 _ TT^" 

36 " y 

In our calculations we neglect the effect of the quarks {Np = 0) and solve the Eq. (fT3|) 
numerically. In what follows we present for the first time the predictions for the total 7*7* 
cross section considering this improved NLO kernel. 

The next-to-leading order corrections discussed above are only part of the corrections 
to the total 7*7* cross section. In a full calculation one should also consider the NLO 
corrections to the impact factors. In the last few years, the real and virtual corrections 



which contributes at NLO has been estimated 1261 and recently the first numerical results 

n 

were presented [27|. These preliminary results indicate that the NLO corrections tend to 
decrease the value of the impact factors. In the phenomenological analyzes what follows we 
will assume that the main energy-dependent NLO corrections comes from the NLO BFKL 
kernel rather than from the photon impact factors. Consequently, our estimates for the total 
cross section at NLO should be consider an upper bound. We will use in our NLO BFKL 
calculations the impact factors as given in Eqs. ([2]) and ([3]) and ujiQ"^, u) taken in the NLO. 
Moreover, following [l^ we consider the Yennie gauge, where ^ = 3. 

Let us now introduce the main formulas concerning the two-photon interactions in the 
color dipole picture. At high energies, the scattering process can be seen as a succession 
on time of two factorizable subprocesses: i) the photon fiuctuates in quark-antiquark pairs 
(the dipoles), ii) these color dipoles interact and produce the final state. The corresponding 
cross section is given by 

n f 

• QlQl) = J2J2 J dz,dz2£r^£r2miruZ^,Ql)\^ (18) 

i,j a, 6=1 
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where i,j are indices associated to the polarization states of the virtual photons and are 
the light-cone wavef unctions of the photon. The variable ri defines the relative transverse 
separation of the pair (dipole) and zi (1 — ^i) is the longitudinal momentum fractions of 
the quark (antiquark) of flavours a,b. Similar definitions are valid for r2 and Z2. The X12 
variable is defined by 



^ Qi + Qi + 4mi + 4m, ^^^^ 



The basic blocks are the photon wavefunction, , and the dipole-dipole cross section, add- 
In the dipole formalism, the light-cone wavefunctions "^{{z, r) in the mixed representation 
(r, z) are obtained through two dimensional Fourier transform of the momentum space light- 
cone wavefunctions "^{{z, k) which are directly related to the impact factors $i discussed 
before (see more details, e.g. in Ref. 28j). The normalized light-cone wavefunctions for 
longitudinally (L) and transversely (T) polarized photons are given by: 

\Hir,z,Q^)\^ = ^5^eJ[.^ + (l-.n4^?(^/r)+mJi^o^(.,r), (20) 

\H{r,z,Q^)\^ = ^5^ej4Q^z^(l-.)^iro^(e,r), (21) 

where Cf and rrif stand for the charge and mass of the quark having flavor / and Kq i are 
the McDonald-Bessel functions. As usual, the quantity e is defined as = z{l — z) Q"^ + rn?^. 
The quark mass mj plays a role of a regulator when the regime of ^ is reached. 
Namely, it prevents non-zero argument for the modified Bessel functions KQ^i{er) towards 
^ 0. 

In Ref. [ITi], Timneanu-Kwiecinski-Motyka (TKM) used the saturation model pro- 
posed by Golec-Biernat and Wusthoff to describe ep collisions [l5| and generalized it two- 
photon interactions at high energies. The basic idea is that the dipole-dipole cross section 
crdd(a^i2, ri, r2) has the same functional form as the dipole-proton one and is expressed in 
terms of an effective radius refr, which depends on ri and/or r2. One has that [l7|. 



(X™(a;i2, reff) = ao 



1 — exp 



(22) 



15|- 



where the X12 variable is given by the Eq. f|T9|) and (Tq = |o"o, with o"o the same as in Ref. 
The last relation can be justified in terms of the quark counting rule. This model interpolates 
between the small and large dipole configurations, providing color transparency behavior. 
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Q' = Q' = 3.5Gev' 




FIG. 1: Energy dependence of the total 7*7* cross section considering distinct approaches for the 
QCD dynamics (Qf = Ql = 3.5 GeV^). 

c"dip ~ r^, at r ^ l/<5sat and constant behavior at large dipole separations r > l/Qsat- The 
physical scale which characterizes the transition between the dilute and saturated system is 
denoted saturation scale, Qsat oc x~^, which is energy dependent. In Ref. [itI three different 
scenarios for has been considered, with the dipole-dipole cross section presenting in all 
cases the color transparency property {add for ri ^ or r2 — > 0) and saturation 
{<ydd o"o) for large size dipoles. We quote also Ref. i29| for further analytical studies 
in the role played by different choices for the effective radius. In what follows, we use 
the model I from [17], where r^^ = r\r\/{r\ + r\)^ which is favored by the 7*7* and 
data. Moreover, in order to extend the dipole model to large X12 it is necessary to taken 
into account threshold correction factors which constrain that the cross section vanish when 
X12 — » 1 as a power of 1 — 0:12. As in Ref. [l7|, we multiply the dipole-dipole cross section 



by the factor (1 — X12Y . A comment is in order here. One shortcoming of the GBW model 
is that it does not contain the correct DGLAP limit at large virtualities. Consequently, 
we may expect that its predictions are only valid at small values of the photon virtualities. 
Therefore, in what follows we only consider photon virtualities up to 15 GeV^. 

In what follows we present our results for the 7*7* total hadronic cross section. The 
quark-box contribution is not included in our calculations. Moreover, as the LO impact 
factors are dependent, we also estimate the effect in our predictions of a running coupling 
constant. Following Ref. |l^ we assume in the LO calculation that = 0.208. This same 
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value is used in the impact factors in the NLO BFKL calculations. The factor scaling the 
energy s in the ressumed logarithms, sq, is a free parameter not predicted by the theory. 
Here we assume that Sq = Q^, where = \/ Q\ Q\, and that is the scale of a^, when 
the running is considered. 

Let us proceed in comparing the different predictions for the total 7*7* cross section as a 
function of the variable Y = \n{s^^/Q'^). For instance, for Y about 10-11, the corresponding 
7*7* center of mass energy reaches 1 TeV, whereas for Y = 7 — 8 one has W^^ ~ 500 GeV. As 
a first analysis, in Fig[T]we consider photons sharing equal virtualities, Qf = Q2 = 3.5 GeV^. 
The LO BFKL result (double dashed - dot line) presents a steep increasing as a consequence 
of the large effective intercept, (ajp — 1) = ^"^^'^ ln2 0.5. Clearly, this behavior is 
ruled out by the L3/0PAL experimental data at high energies pLJ]. For the NLO BFKL 
approach, the BLM scheme and KMRS approach are considered. We have tested the effect of 
running coupling constant within the impact factors and verified that the magnitude of the 
effect is similar in both cases. Consequently, in what follows we only present the predictions 
considering the running coupling for the BLM scheme. The fixed as = 0.208 result is labeled 
by dashed line, whereas the running coupling calculation is the double dot-dashed curve. 
The effect reaches 10% in the overall normalization at high energies in the case of equal 
photon virtualities and the energy behavior seems to be unchanged. Moreover, we have 
that the KMRS prediction (dot-dashed line) is similar to the BLM one at this virtuality. 
Finally, the solid curve shows the TKM model, which describe the low energy L3/0PAL 
data (F < 5). 

The saturation model and the NLO BFKL approach have distinct energy behaviors. The 
underlying dynamics could be disentangled whether a sufficiently large energy collider to 
be available. For instance, at W^^ ~ 500 GeV the deviation is a factor 2 and reaches a 
factor 3 at VT^^ ~ 1 TeV. The NLO BFKL has a flatter energy growth in contrast with the 
saturation model. Roughly, at high energies one has a oc (s^^/Qi(52)°'^^ for the TKM model 
and a oc {s^^/QiQ2)^''^ for NLO BFKL with fixed coupling constant in impact factor. An 
important remark is that although the energy growth of the total cross section predicted by 
the BFKL equation is slower when the NLO corrections are included, its solution still present 
the diffusion property, i.e. the random walk in the transverse momentum. As the mean width 
of this random walk increases as Vln s, we have that at large s the solution eventually enters 
into the non-perturbative region. Thus, independently how large are the transverse scales 
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= 14.0 Gev', q/ = 3.5 Gev' q/ = = 14.0 Gev' 




(a) (b) 

FIG. 2: Energy dependence of the total 7*7* cross section considering different approaches for the 
QCD dynamics, (a) Ql = 3.5, Ql = 14.0 GeV^ and (b) Ql = Ql = 14.0 GeV^. 

of the scattering objects, there is always an energy beyond which perturbative theory loses 
its predictive power. 

In Fig. [2], different configurations for photon virtualities are taken into account: (a) the 
assymetric case Q\ = 3.5 and = 14.0 GeV^ and (b) the symmetric case Qi=Q2 = 14-0 
GeV^ with a somewhat larger virtuality. In both cases the deviations are stronger than the 
previous analysis. At F ~ 11 the saturation model can reach a factor 6 above the NLO 
BFKL (BLM) in plot (a) and a factor 7 in plot (b). Notice the small effect of running 
coupling constant in NLO BFKL when at last one of virtualities is large. The effective 
power for the saturation model has increased in such way that a oc {s^y/QiQ2)^''^'^ and for 
the NLO BFKL (BLM) it remains unchanged. Therefore, the difference between the BFKL 
NLO (BLM) and TKM predictions grows with Q^. In contrast, the virtuality dependence 
obtained using the KMRS approach is distinct, which implies a steeper energy dependence 
for the NLO BFKL (KMRS) prediction. In this case, we have that at F ~ 11 the saturation 
model can reach a factor 2.5 above the NLO BFKL (KMRS) in plot (a) and a factor 3 in 
plot (b) . Concerning the saturation model, we have also checked that the main contribution 
comes from the small size dipole configurations where add oc r^g and saturation becomes 
incleasingly important at W-y-y ~ 1 TeV. Therefore, these effects must can be strong for the 
photon structure funtion, F2^(x, Q"^), where one of photons is real. 
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As a conclusion, the results presented here will be important at possible future e"'"e~ 
colliders (TESLA,CLIC), where the available energies might reach 500 GeV or even 1 TeV. 
In that regime the underlying dynamics could be disentangled. On the other hand, it is 
interesting the comparison of the present results with those ones using BFKL approach 
corrected by sub-leading effects p^l] (where the energy behavior seems to be similar to TKM 



model) and Regge phenomenology for 7*7* interactions as the two-Pomeron models 19|. 



Acknowledgments 

VPG thank G. Chachamis and A. Sabio Vera for useful discussions. This work was 
partially financed by the Brazilian funding agencies CNPq and FAPERGS. 



[1] A. H. Mueller, Nucl. Phys. B415, 373 (1994); A. H. Mueller and B. Patel, Nucl. Phys. B425, 
471 (1994); N. N. Nikolaev and B. G. Zakharov, Z. Phys. C49, 607 (1991); Z. Phys. C53, 331 
(1992) . 

[2] L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); E. A. Kuraev, L. N. Lipatov, V. S. Fadin, 
JETP 45, 1999 (1977); I. I. Balitskii, L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 

[3] I. Balitsky Nucl. Phys. B 463, 99 (1996); Y. V. Kovchegov, Phys. Rev. D 60, 034008 (1999); 
Phys. Rev. D 61, 074018 (2000). L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 
2233 (1994); E. lancu, A. Leonidov, L. McLerran, Nucl. Phys. A 692, 583 (2001); E. Ferreiro, 
E. lancu, A. Leonidov, L. McLerran, Nucl. Phys. A 703, 489 (2002); J. Jalilian-Marian, A. 
Kovner, L. McLerran and H. Weigert, Phys. Rev. D 55, 5414 (1997); J. Jalilian-Marian, A. 
Kovner and H. Weigert, Phys. Rev. D 59, 014014 (1999), ibid. 59, 014015 (1999), ibid. 59 
034007 (1999); A. Kovner, J. Guilherme Milhano and H. Weigert, Phys. Rev. D 62, 114005 
(2000); H. Weigert, Nucl. Phys. A703, 823 (2002). 

[4] G. P. Salam, Nucl. Phys. B 461, 512 (1996); A. H. Mueller and G. P. Salam, Nucl. Phys. B 
475, 293 (1996). 

[5] A. H. Mueller, Phys. Lett. B 396, 251 (1997). 

[6] V. S. Fadin and L. N. Lipatov, Phys. Lett. B 429, 127 (1998); G. Camici and M. Ciafaloni, 
Phys. Lett. B 412, 396 (1997) [Erratum-ibid. B 417, 390 (1998)] 



12 



[7] G. P. Salam, Acta Phys. Pol. B 30, 3679 (1999). 

[8] S. J. Brodsky, V. S. Fadin, V. T. Kim, L. N. Lipatov, G. B. Pivovarov, JETP Lett. 70, 155 
(1999). 

[9] S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983). 
[10] Y. V. Kovchegov and A. H. Mueller, Phys. Lett. B 439, 428 (1998). 
[11] R. Nisius, Phys. Rept. 332, 165 (2000). 

[12] J. Bartels, A. De Roeck, H. Lotter, Phys. Lett. B 389, 742 (1996); J. Bartels, C. Ewerz, 
R. Staritzbichler, Phys. Lett. B 492, 56 (2000); A. Bialas, W. Czyz, W. Florkowski, Eur. 
Phys. J. C 2, 683 (1998); J. Kwiecinski, L. Motyka, Eur. Phys. J. C 18, 343 (2000); N. N. 
Nikolaev, B. G. Zakharov, V. R. Zoller, JETP 93, 957 (2001); S. J. Brodsky, F. Hautmann, 
D. E. Soper, Phys. Rev. D 56, 6957 (1997); Phys. Rev. Lett. 78, 803 (1997). 



[13: 

[14 

[15: 

[16 



[17: 

[18 

[19 
[20 

[21 

[22' 
[23 



M. Boonekamp, A. De Roeck, C. Royon, S. Wallon, Nuc. Phys. B 555, 540 (1999). 

S.J. Brodsky, V.S. Fadin, V.T. Kim, L.N. Lipatov, G.B. Pivovarov, Pis'ma ZHETF 76, 306 

(2002) [JETP Letters 76, 249 (2002)]. 

K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 60, 114023 (1999); Phys. Rev. D 59, 014017 
(1998). 

J. Bartels, K. Golec-Biernat and H. Kowalski, Phys. Rev D66 (2002) 014001; H. Kowalski 
and D. Teaney, Phys. Rev. D 68, 114005 (2003); E. lancu, K. Itakura and S. Munier, Phys. 
Lett. B 590, 199 (2004). 

N. Timneanu, J. Kwiecinski, L. Motyka, Eur. Phys. J. C 23, 513 (2002). 

N. N. Nikolaev, J. Speth and V. R. Zoller, Eur. Phys. J. C 22, 637 (2002); J. Exp. Theor. 

Phys. 93, 957 (2001) [Zh. Eksp. Teor. Fiz. 93, 1104 (2001)]. 

A. Donnachie, H. G. Dosch and M. Rueter, Phys. Rev. D 59, 074011 (1999). 



E. Gardi, J. Kuokkanen, K. Rummukainen and H. Weigert, arXiv:hep-ph/0609087 



Y. V. Kovchegov and H. Weigert, [aFXiv:hep- ph/0609090; I. Balitsky, [arXiv:hep-p h/0609105l, 
G. Abbiendi et al. [OPAL Collaboration], Eur. Phys. J. C14, 199 (2000); M. Acciarri et al. 
[L3 Collaboration], Phys. Lett. B519, 33 (2001). 
L. N. Lipatov, Phys. Rept. 286, 131 (1997). 

D. A. Ross, Phys. Lett. B 431, 161 (1998); M. Ciafaloni, D. Colferai and G. P. Salam, Phys. 
Rev. D 60, 114036 (1999); E. Levin, Nucl. Phys. B 545, 481 (1999); M. Ciafaloni, D. Colferai, 
G. P. Salam and A. M. Stasto, Phys. Rev. D 66, 054014 (2002). 



13 



[24] G. P. Salam, JHEP 9807, 019 (1998); M. Ciafaloni, D. Colferai and G. P. Salam, Phys. Rev. 
D 60, 114036 (1999). 

[25] V. A. Khoze, A. D. Martin, M. G. Ryskin and W. J. Stirling, Phys. Rev. D 70, 074013 (2004) 
[26] J. Bartels, S. Gieseke and C. F. Qiao, Phys. Rev. D 63, 056014 (2001) [Erratum-ibid. D 

65, 079902 (2002)]; J. Bartels, S. Gieseke and A. Kyrieleis, Phys. Rev. D 65 (2002) 014006; 

J. Bartels, D. Colferai, S. Gieseke and A. Kyrieleis, Phys. Rev. D 66, 094017 (2002). 
[27] J. Bartels and A. Kyrieleis, Phys. Rev. D 70, 114003 (2004); G. Chachamis and J. Bartels, 

Proc. of the Int. Workshop on Diffraction in High-Energy Physics (DIFFRACTION 2006), 

PoS (DIFF2006) 026. 

[28] V. Barone and E. Predazzi, High-Energy Particle Diffraction, Springer- Verlag, Berlin Heidel- 
berg, (2002). 

[29] C. Marquet and R. Peschanski, Phys. Lett. B 587, 201 (2004). 

[30] J. Kwiecinski, L. Motyka, Phys. Lett. B 438, 203 (1998); Acta Phys. Pol. B 30, 1817 (1999). 



14 



